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1 Results 



In this note, we consider the free Bosons hopping on networks (general graphs) 
and we consider the condition for Bose-Einstein condensation. Let T = {V,£} 
be an infinite graph where V is the set of vertices and £ is the set of edges. 
Here for simplicity, we assume that a pair of vertices i and j is not connected 
by more than two edges (no multiple edge), thus we can denote the unique edge 
connecting vertices i and j by (i, j). The graphs we consider is connected and 
not oriented. Furthermore for simplicity of presentation, we assume that the 
graph does not contain any self-loop, i.e. (i, i) ^ £. 

Next we introduce free Bosons hopping on vertices of T. Let dj and a* k ( 
j, k in V) be the creation and annihilation operators satisfying the canonical 
commutation relations, [aj,a£] = 5jk- Smeared Boson operators, a*(f) and 
a(/), are defined as 

«*(/) = £ a *y* , «(/) = E a *7 fc (i-i) 

jev jev 

where fj is a complex function on V. The function fj is referred to as a test 
function. Obviously, 

[a(f), a*{g)] = {f,g)i2 {v) l 
where (f,g) is the inner product of l 2 (V), 

if,9h(V) = 

jev 

To be definite, we choose the set of rapidly decreasing test functions on V as 
the test function space. Rapidly decreasing functions on V can be introduced 
as follows. By definition, the graph distance dist(j, k) of two vertices j, k is the 
minimum of the number of edges giving rise to a path connecting j, k. Fix the 
origin o of the graph and a complex function / on V is called rapidly decreasing 
if 

^rfisi(i,o) n |/, | < oo 

-jev 

for any positive integer n. The set of all rapidly decreasing functions is denoted 
by S. Set 

W(f) = eM^(a*(f) + a(f))). 

Then, 

W { f)W {9 ) = eM-i Im{{f t 2{V)) )W(f + g), 
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where Im({f,g)p)(y) is the imaginary part of the inner product of l 2 (V). We 
denote 2l(r) by the Weyl CCR C*-algebra generated by W(f) (/ e S). 

Now we consider the free Hamiltonians on our network. On Euclidian spaces, 
the time evolution of the free Boson is determined by the Laplcian acting on 
the test function space. Two natural candidates of the graph analogue of the 
Laplacian are often used. One is the discrete Laplacian and the other is the ad- 
jacency matrix. By degree of a vertex j we mean the number of edges connected 
to j and we denote the degree of j by d(j). (The degree is sometimes called 
coordination number.) The discrete Laplacian — Ar on T is defined by 

-A r f(j) = d(j)f(j)- /W' (1-2) 

k:(k,j)d£ 

where the sum is taken for all vertices k connected to j by an edge of T. The 
discrete Laplacian is a positive operator on 1 2 {V). The adjacency matrix At is 
the off-diagonal part of the discrete Laplacian, 

Arf(j)= J2 /(*)■ (1-3) 

k:(k,j)e£ 

Thus when the degree of the graph is a constant d(o), —A T + d(o) = —A and 
there is no physical difference of these two operators. When the degree is not 
constant, we regard the adjacency matrix — A T as a Schrodinger operator on a 
graph and the degree plays a role of potential. 

We introduce assumptions for graphs to describe our main result. 

Assumption 1.1 (i) We suppose that the Graph V is connected and that T is 
the limit of an increasing sequence of finite connected subgraphs T n = {V n ,£n}, 

r = u~ =1 r„ = lim r„. 

n^oo 

(ii) We assume the F0lner condition is valid: 

l im = . (1.4) 



n— »oo 



where \T n \ is the number of vertices ofT n and \dT n \ is the size of the boundary 
ofT n . More precisely, \dT n \ is the number of vertices ofT n connected to the 
complement ofT n by an edge ofT. 
(in) The degree of the graph is bounded, sup<i(j) < oo. 

The conditions (i) and (ii) above are referred to as van Hove condition in 
statistical mechanics and they are used to handle thermodynamical quantities of 
infinite volume systems. The condition (iii) ensures boundedness of the discrete 
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Laplacians and adjacency matrices. We note that the condition (ii) excludes cer- 
tain non-amenable graphs such as Cayley trees while the condition (iii) excludes 
the complete graph (all pairs of vertices connected by edges) corresponding to 
the mean field free model. The case of Cayley trees is studied in [T3|. 

For a finite sub-graph r„, the discrete Laplacian and the adjacency matrix 
are denoted by — Ar„ and Ar n - If the graph is finite, the constant function 1 
is the unique ground state for the discrete Laplacian, 

- A r 1 = 0. 

The unicity (up to multiplicative constant ) of the ground state follows from 
the Perron-Frobenius theorem for positive matrices. 

The quantum Hamiltonian is the second quantization dT of these operators, 
thus H = dT(— Ar) or H = aT(— Ar). In what follows, we consider mainly the 
second quantization of the discrete Laplacian and we set 

H r = dT(-Ar) = - \ ]T (a*a k + ofo-)}. (1.5) 

jev k-.(k,j)e£ 

Then, 

[Hr,a*(f)] = a*(-A r f) 

The time evolution of observables in the Weyl CCR C*-algebra -4(r) on T 
is determined by 

a T t {W{f)) = W{e-^f). 

For each finite sub-graph T„, let De the equilibrium state(grand canonical 

ensemble) at the inverse temperature /3 with the mean density p: 

^ ) W))=exp(-i(/,i±^/ )) 

V<ff\a*(f)a(g)) = (g, *_?f** f) (1-6) 

where z n (0 < z n < 1) is determined by 

z n is the exponential of the chemical potential \x n , z n = exp(— f3/j, n ). 

When T is a finite graph, we define Zg(V) = {(fj) G 1 2 {V) \ J2 je v h = 1- 
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By E we denote the projection from l 2 (V) to Iq(V). We set tr (Q) = tr(E QE ) 
and 

1 zeP Ar " 

p{P)= sup limsup—r tr(- g^r)- (1-8) 

0<z<l n |r n | 1-Ze^ r ™ 



Then, 



because 



1 e' 3Ar " 
p(/3) = limsup— tr (- ^— ) (1.9) 



ze /3A r „ 

tr (^ — __ fl A. ) = i — I + tr o 



The following can be shown in the same way of the Bose gas on the Euclidean 
space. ( See |2 dU and [7]. ) 

Proposition 1.2 (i) Suppose that~p((3) is infinite. Then, for any p > 0, there 
exists Zoo and a subsequence of finite graphs T n u\ such that lirm z n u-\ = Zoo < 1 
and 



(p.p), .. <» n\ , 1,. 1 + z 00 e' 3Ar 



\w(f)) = lim^(W(/)) = exp(--(/, l°° /)). (1.10) 



^VCfMs)) =lim^/ ) (a*(/)a( 5 )) = (g, - ^^ /). (l.H) 

n l — z oc e p " r 

fnj Suppose that~p(/3) is finite. 

(iia) If p < ~p(f3), there exists lim„ z n = Zoo < 1 and a subsequence of finite 
graphs T n u\ such that the equations Hl.lUp and HI. lip are valid. 

(iib) Suppose p > ~p((3). Take any subsequence of finite graphs T n ^ satisfying 

1 e^r-w 
p((3) = lim 1 tr gs^— 1.12 

Then, 

limz n(l) = 1 , lim— -=p-p{f3). (1.13) 

Furthermore, suppose that any f = (/j) wii/i compact support is in the domain 
of {— Ar) -1 / 2 and that the following uniformly boundedness of matrix elements 
holds. 

sup (ft, . . ft) < oo. (1.14) 
i,j'ev l-^rj 
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Then, any rapidly decreasing f and g are in the domain of (— Ar) (hence 
in the domain of (1 — e /3Ar ) _1 / 2 y ) and, 



where 



^<>\w{f)) = \im^f{W{f)) 



exp(--(p - p(/3))|xr(/)|^) exp(—(/, ^^ Z)) (1.15) 



^' p) (o*(/)o(ff)) = ]im^f\a*(f)a(g)) 



e /3A r 



(P - p(/?))xr(/)Xr(5) + (.9, 1 _ g/3Ar /) (1-16) 



Xr(/) = E fi- ( L17 ) 



We say Bose-Einstein condensation occurs for the second case, p > p(f3), 
in the above proposition. Due to the equation (|1.15|) and (|1.16l) . we have off- 
diagonal long range order and the U(l) gauge symmetry breaking for the state 
<Pj? with a high mean density in the following sense: 

lim ^ P \a* aj ) = ip-m) * = 4 M (a*)4 M ( aj ). 

di8t\%,j) — ►co 

(pj? ,ri ^ is decomposed into the following factor states. 

4 M (W(f)) = eM^-p{ti)) 1/2 ^ e Xr{f) + e-^Xv(f))} 

xexp(- i (/, T ±-^-/)) (1.18) 

In another word, >p^' p ^ 18 an integral of ip^?' 11 ' 6 ^ as functions of rj and 8. 

When the one-particle Hamiltonian is a discrete Laplacian, we present the 
condition for Bose-Einstein condensation in terms of the simple random walk on 
r. By the simple random walk, we mean that the random walk can move only 
to adjacent vertices and that the probability of the jump from one vertex i to 
another adjacent vertex is l/d(i). The probability of the random walk moving 
to adjacent vertices are all equal. 
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Definition 1.3 Letpisr(j) be the probability of the random walk starting from j 
and returning to j at the Nth step for the first time. Let 9jv(j) be the probability 
of the random walk starting from j and returning to j at the Nth step. 
The simple random walk on V is recurrent if ^2^ =2 PN(j) = 1 for each vertex 
j. The simple random walk on T is transient */X)w=2^(i) * s finite f or each 
vertex j . 

It is well known that the random walk is not recurrent if and only if it is 
transient, (c.f. |14p For general inhomogeneous networks, we use stronger 
conditions to describe our results on Bose-Einstein condensation. 

Definition 1.4 (i) The simple random walk on T is uniformly transient if 

oo 

sup Y] q N (j) < oo. (1.19) 

(ii) The simple random walk on T is uniformly recurrent if for any large K there 
exists M such that for any j 

M 

K <J2 8* C?)- (1-20) 

JV=2 

We do not know whether uniformity conditions of p. 19(1 and 11.20JI are used in 
other contexts of probability theory. The same conditions may be expressed in 
other terms. 

Theorem 1.5 We consider the free Bosons on a graph T satisfying previous 
conditions Assumption 1.1 and assume that T contains no self-loop and no mul- 
tiple edge. Let the one-particle Hamiltonian be the discrete Laplacian — Ar- 

(i) Suppose that the simple random walk on T is uniformly transient. Then, 

p(0) < oo 

and any rapidly decreasing f is in the domain of (— A) -1 / 2 , hence, alll the 
assumptions of Proposition^^ (iia) and (iib) are valid and the Bose-Einstein 
condensation occurs. 

(ii) Suppose that the simple random walk on T is uniformly recurrent. Then, 

p([3) = oo, 

and the conclusion of Proposition^^ (i) holds. 

Example 1.6 (Periodic lattice) We consider an infinite graph T on which 
the abelian group Z" acts graph automorphisms. By we denote this Z" ac- 
tion. We assume that 
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(i) the isotropy group is trivial , namely, for any vertex i, Tk(i) = i implies 
fc = 0. 

and 

(ii) the Z" action is co-finite in the sense that the quotient graph Tq = T/Z v is 
a finite graph. 

A graph satisfying these conditions (i) and (ii) is referred to as a periodic lat- 
tice and To will be called the fundamental domain ofT. In periodic lattices, by 
definition, the simple random walk is uniformly transient (resp. uniformly re- 
current) if it is transient (resp. recurrent). It is known that the random walk is 
transient if and only if the dimension v is greater than or equal to three, 3 < v . 
See \10f . In the periodic case, we can replace lim sup with lim in taking infinite 
volume limit and we do not have to take subsequences of finite graphs. 

Example 1.7 (Defects) Now we remove edges from graphs and we call re- 
moved edges defects. We denote the set of defects by T>. We assume that the 
density of defects vanishes 

llm iv I = °- 

Then, our proof shows that if the simple random walk is uniformly transient for 
the initial graph, BEC occurs for the graph with defects. If the density of defects 
is positive, further consideration is required and we do not know condition for 
occurrence or absence of BEC in this case. 



Next we discuss the case when the one-particle Hamiltonian is an adjacency 
matrix. The transiency of the random walk does not determine the occurence 
of BEC. Such examples are provied in 0] and|S]- The examples discussed in 0] 
and|S] are star graphs and comb graphs for which the simple random walk are 
recurrent . (See In 0], I- Brunelli, G. Giusiano, F. P. Mancini, P. Sodano 

and A. Trombettoni considered the star graph and the adjacency matrix has 
a spectral gap. In such situation, it is easy to see the existence of the critical 
density and occurrence of BEC. In 0, R.Burioni e D. Cassi, M. Rasetti, P. 
Sodano and A. Vezzani investigated BEC on the comb graph. The vertices of 
a comb graph are same as Z 2 but some edges are removed from Z 2 . In our 
terminology, the density of defects is strictly positive. 

When the graph is a periodic lattice and the adjacency matrix is the one- 
particle Hamiltonian, we have the same results as the integer lattice, so BEC 
occurs if and only if 3 < v. For the continuous space case on R^, this follows 
from a result of W.Kirsch and B.Simon in [S] (See also • W.Kirsch and 
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B.Simon considered the Schrodinger operator with a periodic potential and they 
proved that behavior of the density of states in the vicinity of the bottom of 
energy is same as that for free Laplacian. The argument of W.Kirsch and 
B.Simon works for our periodic lattice as well. We obtain the same result when 
we add a periodic potential to the one-particle Hamiltonian. 
We set 

^n=^\k t \<n (i=l,2,...i/)Tfc(r ) (1-21) 

(p) 

By T„ , we denote the graph obtained by the periodic boundary condition on 
r! p ' . The vertex set of is same as that of T n and the edge set of ^ is the 
union of T n and additional edges in the way that there is a natural (Z2 n +i) l/ 
action on compatible with the shift of T. More precisely, if i £ V n , j ^ V„, 
£ £ and there exists m £ Z" such that \m\ = 1 and r m (i) — k for some k 
in V n , (i, k) is an edge of T„\ Using this periodic graph rl P \ we can introduce 
the one-particle Hamiltonian with the periodic boundary condition. 

Theorem 1.8 Suppose the graph T is a periodic lattice and v(j) is a (real) 
periodic function on the set of vertices V. Let E n (resp. E) be the supremum of 
the spectrum of A r ( P ) —v (resp. At — v) acting on l 2 (V n ) (resp.l 2 (V) ) . Consider 

h„ = E n - A r(p) + v, (1.22) 



h = E-A r + v. (1.23) 

Let ip { r?' p) be the equilibrium state associated with the second quatization of h n 
with a mean particle density p. It is determined by 

e p) (^(/))-exp(-l (/ , | + ^:^ /)), 
±-<pW*){N n )=p. (1.24) 

(i) Suppose that v < 2. BEC does not occur. For any p the limits lim„ z n = 
z < 1 and 

1 1 -I- 7p~P h 

lim ^Hw(f)) = V {M {W{f)) = exp(--(/, -— /)) 

n—i-oo 4 1 — ze p 

exist. (p(l 3 -'P'> is a translationally invariant factor state o/2l(r). 

(ii) Suppose that v > 3. Then the rapidly decreasing functions are in the domain 
of h~ x / 2 and we can define the critical density pp via the trace tr ; 2( Vo ) over the 
fundamental domain of our periodic lattice: 

1 e~ (3h 1 ^ e _/3/l 

aeVa 
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Suppose pp < p. Then, lim n _, oc> z n = 1 and 

^ p \w(f)) = lim^f\w(f)) 

1 11+ (=-0 h 

= exp(--(p - p((3))\xAf)\ 2 ) «P(- Z (/, 73^/)) ( L26 ) 



< ,w (o*(/)a( fl )) = (p- p(J3))Xr{f)Xr(9) + Cff> 1 _ e -p h f) 



where 



Xr(/)= Hm |r„| 1/2 (0„,/) i2(Vii) . (1.27) 

n — >oo 

and f2„ is the positive normalized ground state of h n uniquely determined by 

h n VL n = 0, (n n ,O n )(2(v„) = 1. (1.28) 

Alternatively, Xr(f) «s described by the inner product of the positive periodic 
ground state fl — {fl(i)} for the infinite volume Hamiltonian h. Let Q belong to 
l°° satisfying the normalization condition 

El^)! 2 = 1 (1-29) 

ieVo 

Then, 

*r(/) = = (nj)mv), (1-30) 

which converges for any rapidly decreasing f. 



Next we explain key points in our proof of Theorem 11.51 and 11.81 The first 
point is finiteness of the critical density 7j(/3). For x > and < z < 1, 

^° _ \ m — mfjx 

1 - ze~P x ~ ^ 

m— 1 

Hence, to examine finiteness of ~p((3), we consider the trace of Gibbs weight, 
jjT-j- tr(e - ^ "•). Let Pp(v n ) be the projection l 2 (V) from to l 2 (V n ) and d n (resp. 
d) be the degree of T„(resp. T). By Trotter-Kato product formula we will see 

ze -f3h„ ze -f3h n 
toi'(v n )( ,_ zp - hJ ^ tr mv)(Pi*(r n ) 1 _ phn Pp(r.)) (1-31) 
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Then, by use of the F0lner condition ljl.4|l . we show that we can replace d n with 
d in the right-hand side in the above inequality. Then, we obtain 

Pi/3) < sup (Su Z€ k] < C U, -—± r8i) (1.32) 

^ev\ l-ze^ J ;2(v) V h+{l-z) J , 2(v) 

for a positive constant C. If the diagonal of the Green function is bounded 
uniformly in i, we obtain finiteness of the critical density ~p(f3). It is east to see 
that boundedness of the diagonal of Green function is equivalent to transiency 
of the simple random walk. This shows the first part of Theorem II .51 Another 
inequality can be obtained by similar way and the divergence of the Green 
function on the diagonal is implied by recurrence of the simple random walk. 

When the graph is a periodic lattice, the situation is simpler as we employed 
the periodic boundary condition. It is easy to see 

1 zer^ hn 

„'iS,]rJ tr "' r -» ( i-.e-». ) 

Thus, proof of finiteness or divergence of the critical density reduces to finiteness 
of the diagonal of Green functions again. We show an inequality between the 
Dirichlct forms attached to our discrete Schrodinger operators, (c.f. Lemma 
I3.9fl . This tells us that the behavior of the density of states at the bottom 
of spectrum for our discrete Schrodinger operators with a periodic potential is 
same as that for our discrete Laplacian. 

In Theorem 11.51 we assumed a number of technical conditions on graphs. 
Conditions of as no self loop and of no multiple edge are not essential. If 
the graph has self loops we have only to consider random walks which stay 
on the vertex with certain probability. To handle multiple edges, we modify 
the definition of the adjacency matrix and its entry is the number of edges 
connecting vertices labeling entries. On the other hand, the F0lner condition 
(|1.4f) plays an crucial role in our argument. Even though we are unable to 
prove Bose-Einstein condensation without using the F0lner condition, we believe 
the same result is valid for graphs without our F0lner condition. One such 
example is the Cayley tree associated with the free group. In ^3] M.van den 
Berg, T.C.Dorlas and V.B.Priezzhev proved Bose-Einstein condensation for the 
free Boson hopping on a Cayley tree. Their argument is based on explicit 
computation of the spectrum of the discrete Laplacian. 

In Section 2 we present our proof of Theorem 11.51 and Section 3 is devoted 
to Theorem Ol 
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2 Proof of Theorem II. 5L 

In this section, we prove Theorem 11.51 Let W be a subset of vertices V Let 
us regard l 2 (W) as a closed subspace of l 2 (V) and we denote the orthogonal 
projection from l 2 (W) to l 2 (V) Pw- We regard the degree d(j) of F as a multi- 
plication operator on l 2 (V) denoted by d, while the degree of the subgraph r„ 
is denoted by d n . d n is different from d{j) on the boundary of T n . We set 

d = supiirO) < co. 
jeV 

By trs)(Q) we denote the trace of an operator Q on a Hilbert space Sj. For 
simplicity, Ar(resp. Ar n is denoted by A (resp. A n ). 



Lemma 2.1 

tr 



( v„)(e^ r ") 

< tr l 2 {v) (P dTn f3Ae^P er J + tr !2(Vn) (exp(P rf n t /3 A r P r .„»0) 

< tr lHV) (P aTn pAe^P d r n ) +tv lHv) (e^ P dr J 
+tr /2(v) (P v „e' 3Ar P v J (2.1) 

where T™* is the interior ofT n , T" lt = V n — dT n . 

Proof. The positive discrete Laplcian — Ap„ is composed of the diagonal term 
(multiplication operator) d n and the off-diagonal part identical to the adjacency 
matrix A n of. The difference of the (positive) discrete Laplcian — Ar„ and 
— Pr„ Ar Pr„ is the diagonal part. Using the Trotter-Kato product formula, we 
obtain 



tr i 2 (Vji) (e" /5( " Ar " ) ) = lim tr, 



N — >oo 



(1 



pA n , 
N ■ 



-0 l - 



(2.2) 



Set 



X(N) 



(1 



N 



X(N) = P rint X(N)P T i 



Then, the right-hand side of the equation l|2.2|l is 

tr l2{Vn) (X(N) N ) = tx lHVn) (X(N) N P d r n ) +tx PiVn) (X(N) N P r 
tT P( y n) (X(N) N P aTn ) + trpty^XiN^XiN)) 
+tT l 2 iVn) (X(N) N - 1 P r ^ t X(N)P dr J 
= tr lHVn) (X(N) N P drn ) + trp^XiN^XiN)) 

+tr(X(NY 



N 



(2.3) 



12 



As all the terms in (|2.3(l are the trace of non-negative matrices, 

tr P(Vn) (X(N) N ) 

< tv PiVn) (X(N) N P dTn ) +tr P{Vn) (X(N) N - 1 X(N)) 

+ ^tr P{Vn) (X(N) N - 1 P r ^pAP 9rn e^ d -) 

< tr lHVn) (X(N) N P dT J+tr lHVn) (X(N) N - 2 X(N) 2 ) 
+ ^tr P{Vn) (X(N) N - 1 P rM pAP drn e-^) 

< ti l 2 (Vn) (X(N) N P drn )+tr l 2 (Vn) (X(N) N ) 
+tr l 2 {Vn) (X(N) N - 1 P rTt [3APgr n ). (2.4) 

Now we take N to infinity in (|2.4|) , and we obtain 

tr, a(VB) (e^ r ») 

< tr P( y n) (e^P d r n ) + ^^{e^^^) 
+tY P{Vn) {e^P T ^l3AP aTn ) 

< tr P[Vn) {e^P dTn ) + tr lHTtt) (e 0P ^ tA ^ P ^) 
+tY P(Vn) (e^pAP dTn ). (2.5) 

Next we show that for any subsets of vertices Wi and W2 satisfying W\ C W2 C 
V, 

tr i2(VVl) (e^ Ar » p -i) < tr P(W2) (e^ A -" p ^) < tr P(Wa) (e^ Ar ) (2.6) 

Again, this is due to the Trotter-Kato formula and non-negativity of matrix 
elements of the adjacency matrix. 

tr P(Wl) (e^ A r„P Wl) = lim tr w ((Pwi( i + ^j^-M./")*) 

iV — >-oo IV 

< lim trp (W2) ((P W2 (l + ^)P W2e ^"^) w )=tr i2(W2) (e^ A -^ 2) 

< lim tr i2(W2) (((l + ^)e~^f) = tr ;2(W2) (e~ 

IV — >oo IV 

As n is arbitrary, we take n to oo and we obtain, 

tr ;2(W2) (e^ A ^) < tr /2(W2) ( e -^ d -^)) = tr ;2(w . 2) (e^) 
Combined with (|2.6|l , the equation 12. 5|) implies 12.1(1 . End of Proof. 
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Lemma 2.2 There exists a constant C such that 

t ^(Vn)( I j^~ ) < C\dr n \(l3^^ + I ^- z ) + tr lHV) (P rn i g^r PrJ 

(2.7) 



Proof. As the heat kernels are contractive, e^ Ar and e^ Ar ™ have the norm 1. 
By definition, 

tr P(v) (ParJ = |0r n |. 

Thus, we have 

|tr P(v) (P arn /Me^ Av P ar j| < f3\dV n \[3\\A\\ , (2.8) 

and 

tr / 2 (r) (e /3Ar "P ar J < |flT n |. (2.9) 
On the other hand, recall that 



ze 



/3A r 



1 _ ^ e /3A r „ 

M=l 



E 



,M M/3Ar„ 



for < z < 1 and /? > 0. We combine the trace of this geometric series with 
the inequality J2HJ EHl and to get End of Proof. 



Lemma 2.3 (i) Suppose that the simple random walk is transient. Then, the 
delta function 8j is in the domain oJ-Ap 1 ' 2 

(ii) Suppose further that the simple random walk is uniformly transient. Then, 



*<.-W) < Cr if ( 2 - 10 ) 

where rj°' is the probability of the simple random walk starting from i arriving 
at j. 

Proof, (i) Note that 

qN(j)=(s i ,(~Ar) N 8 i ) (2.11) 
V a / p(r) 
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Then, for positive z, we have the following Neuman expansion: 
Si, — -7— Si J = f Si, 



z-A r / ;2(r) V 'd + z-Ar / ;2(r) 

^ V ^ + z d + z y P(r) 

1 

<(l+^ gjv (i) ) <00 . (2.12) 

(ii) Let ry (iV) be the probability of the random walk starting from i and arriving 
at j at the Nth step and let rff (N) be he probability of the random walk starting 
from i and arriving at j at the Nth step for the first time. Then, 

N-2 

r ij (N)=r^(N)+ £ nj (k) Pj (N - k). (2.13) 

k—dist(i,j) 

Set 

OO OO 

nM= E njW* n > ^\z) = E rf(iV)A 

N=dist(i,j) N=dist(i,j) 
00 00 

H z ) = E pjW^. U z ) = E ^w- 

JV=2 JV=0 

Multiplying z N and adding in N we have 

f ij (0) = rg ) (^)+r ij (^.(^) (2-14) 
If \z\ < 1 , fij(z), r\®\z) andpj(z) converge absolutely. By definition, 

and the random walk starting from j is transient if and only if p < 1. As 
we assumed that the random walk is uniformly transient, there exists a positive 
e such that , for any j,pj(l) < 1 — e. As a consequence, 

r (o) (0) 

Returning to the Neumann expansion of the resolvent, we have the following 
relation: 

r N (ij)=(Si,(±A r ) N 6 j ) . (2.16) 
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If i and j are different, we obtain 

^) P(r rl ( < " ( >"^) P(r) ^<f- <™> 



End of Proof. 



Lemma 2.4 Set 

1 

|f ^tr ;2(v) (P v „ T — 
Assume that the simple random walk is uniformly transient. Then, p{(3) is finite. 



p(fi) = limsup — tr ;2(v) (P v „- — PyJ- (2.18) 



Proof. As — Ar is a positive bounded operator, the functional calculus implies 
the following inequality with suitably large C. 



e 0A T x 
<C- 



1 -eP A ? - -A r ' 

By previous lemma, we have 

< |r„| sup^(l). (2.19) 
This implies the finiteness of p(/3). End of Proof. 



Lemma 2.5 Assume that the simple random walk is uniformly transient. Then, 

P(P) < P{P) < oo. (2.20) 

In particular, ~p(f3) is finite. 

Proof. Assuming p([3) < p((3), we show contradiction. Take p satisfying p{[3) < 
p < p((3). There exists a sequence of subgraphs r„( fe ) and Zk such that the 
following is valid: 



z fc e 

lim Zk = < 1. (2.21) 

k 
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However, the inequalities 12.21JI and 1|2.7(1 imply that p < p(f3). End of Proof. 

To prove Theorem 11.51 (ii), we introduce a new notation. For any positive 
integer m we define an augmented boundary: 

d m T n = {ieV„ \dist(i, r£) < m} 

where is the complement of T n in T. Then, obviously, 

\d m T n \ <d rn \dr n \. 

As a result, 

™- > °° Hn| 

Lemma 2.6 Suppose that the simple random walk is uniformly recurrent and fix 
large K and no such that the inequality \1.2(]\) is vaild. Take larger n\ satisfying 

\dnjn\ ^ 1 

rao<ni, | r | < V 2 
/or n > ni . There exists a constant C independent of K such that 

p(P) > CK. (2.22) 



Proof. First by use of Trotter Kato formula, we show 

fee^JOpo;,, > (5i,e^^ p ^5i) lHVn) . (2.23) 
As d„ < P Vn dP Vn and A r „ = Py„^r-Pv„, 



> 



(&,e p r "&)z 2 (v„) = lim (Me «e « ) <5*)z 2 (v„) 

iV — >OD 

lim (<5i,(e e ) &)z 2 (v„) 



Thus we obtained lt2~23|) 
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Next take C\ such that er x — 1 < C\x for x satisfying < x < d. Note that 
e x - z = e x - 1 + (1 - z). By lt2~2"5)l we have 

1 ze /3Ar " 1 z 

jr^'^i-ze/^rj " C W* ri2(V " )( -Ar„+l-^ 
1 

- Cirf^^ 2 (v n )( 



> 



|r„| lKn;v -iV» A r p v „ + i-^' 

i 



1 n| lier n na ni r° 



P v „(-A r )Pv„ 



^(v.) > • (2-24) 



As before we use the Neumann expansion of the resolvent: 
( * < 'iV f .(-Ar)k + l-^ i),a(V " ) 



OO ■[ ! 

^ EW. ^ZTi^^^Ti^^^ 

L — 

* E^s+i^^-^rf+T^^cv-) (2-25) 

as each term in the above Neumann expansion is non-negative. 
As far as i belongs to T n n d ni T^ and L < n% 

= ^jThz^jrh-/ 6 ^ (2 - 26) 

This is because the random walk starting from i cannot reach the augmented 
boundary d ni T n within L steps. 
Now we combine these estimates 

1 ze /3Ar " 
lr^ i2 < v "> ( l-ze^„) 

C\Z ^ l 

- 2ir nd r c l E E d + 1 _ 7 ( Ar d ■ i _ z ) L ^)' 2 (v„) 

(2.27) 
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As a consequence 



§■ inf ft(n(l)) > ^-K. (2.28) 



End of Proof. 

The above lemma completes our proof of Theorem 11.51 
End of Proof of Theorem EU 



19 



3 Proof of TheoremEHl 



In this section, we prove Theorem 11.81 As remarked before, when the one- 
particle Hamiltonian is a Schrodinger operator on Euclidean spaces, the question 
of occurrence of the Bose-Einstein condensation is reduced to behavior of the 
density of states at the bottom of the spectrum. In case of periodic potentials 
W.Kirsch and B.Simon proved that the asymptotic behavior of the density states 
is same as that for the free Schrodinger operator. (See j2] •) The same problem 
is considered in presence of magnetic field by P.Briet, H.D.Cornean, and V.A. 
Zagrebnov in The basic idea of our proof for the periodic lattice is same 
as U], however, there appears some difference between the periodic lattice case 
and the Euclidean case, which we explain below. 

Let us recall that the periodic lattice T = {V,£} is obtained by the funda- 
mental domain Tq = {Vo,£o} and its shift. The shift is denoted by as before. 
The choice of fundamental domain is not unique and we fix one fundamental 
domain Tq — {Vo,£o} satisfying two conditions: 

(i) it is a connected sub-graph of T. 

(ii) the conditions Tk(i) = j, and i,j € Vo imply i = j and k — 0. 
Furthermore, without loss of generality, we may assume 

(iii) r o p) does not possess a multiple edge. 

For our purpose, we may consider a larger block T n as a fundamental domain 
and a smaller group ((2n + l)Zi) v as the shift on our lattice. This is the reason 
why we can assume (iii). 

Now, we identify the vertex set V with Vo x Z" in such a way that the shift 
on Vo x Z" acts via the following formula: 



l 2 (V) is naturally isomorphic to ^ 2 (Vo) <8> l 2 {7i v ), We identify this Hilbert space 
with ; 2 (V ) <gi L 2 {T V ) by use of Fourier transformation F from l 2 (Z") to L 2 (T") 
where the torus T d is identified with [— 7T, 7r]". Our convention for Fourier trans- 
form F and the inner product of L 2 {T V ) are 



Then any operator B commuting with shift is called translationally invariant. 
It is unitarily equivalent to a matrix valued multiplication operator B(p): 



Tk(a,j) = (a,j + k). 





FBF- l f{p) = B(p)f(p) 



f{ P )el 2 {V a )®L 2 {T»). 
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Thus, we have direct integral representation of any translationally invariant 
operator B : 

B= f B(p)dp, 1 2 {V) = f $j p dp 
Jt" Jt v 

where f) p is a |To| dimensional Hilbert space of wave functions satisfying the 
following twisted boundary condition: 

% = {/(a, j) € Z°°(V)|/(a, j + k) = e**/(o, j)}. (3.1) 

(p) 

The above observation is valid for finite periodic graphs r,\ introduce in Section 
1 as well. We have only to replace the 7/ with the finite cyclic group Z„ 
(Z n = Z/nZ). Following the custom of physicists, we call the variable p quasi- 
momentum. 

Ay (p) is a matrix with entires indexed by vertices of the fundamental domain 
Tq. The matrix elements [-<4r]ij(p) of Ar(p) ( for a fixed quasi- momentum p) 
are described as follows. 

Set ei(k,j) — Sij where k € Z" and i,j G r . We regard ej a periodic 
function on V. Obviously ej is in Z°°(V) and {e^z € V )} is a basis of the set 
of periodic functions on V. When the quasi-momentum p is zero, [^4r]ij(0) is 
determined by 

A rej = $^[Ar]tf(0)ei. (3.2) 

As the degree of graph is bounded, the adjacency matrix Ar is a bounded 
operator on l°°(V) and l|3.2fl should be understood as an identity of Z°°(V). 
To consider the case for non vanishing quasi-momentum p, we set 

e?\k,j) = e ik -*6 ij . 

{ef\k, j)} is a basis of the space Sj p of functions with quasi-momentum p. 
Then, the matrix elements of Ar(p)ij is determined by 

Avef = Y,[MiM^ ] - (3-3) 
By construction and translational invariance, we obtain the following formulae. 

Lemma 3.1 (i) If i and j are connected by an edge inside the fundamental 
domain Tq, 

[Ar] ij (p) = [Ar} ij (0)- (3-4) 
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(it) If i and j are connected by an edge bridging adjacent blocks ofT, 

[A r ] ij (p)=e i9 ^[A T } ij (0). (3.5) 

Here Oij{p) is real and 9ij(p) = —Oji(p). It is a linear combination of the 
compnent p(k) of the quasi-momentum p — (p(l),p(2), • • ■ p(v)) described as 
follows: If the edge (i,j) of £q corresponds to an edge (i,j) ofT such that 

i = ((ii,»2i • • • >iu),a) , j = ((ji, • • • ,j v ),b) eV = Z"xV , 

V 

1 < y^Q'n - in) <V, i n < j n (n = 1, 2, • • • , v), 
n=l 

v 

Qijip) = ^2Pn{jn ~ in)- (3.6) 

n=l 

(Hi) If(i,j) is not an edge ofT, 

[A r ]ij(p)=0. (3.7) 



For later convenience, we set 9ij (p) = when i and j are connected by an edge 
inside the fundamental domain To and we can write 

[Ar} ij {p) = e ie ^[Ar}i j (0) (3.8) 

for any i and j. The factor e l6ii ^ corresponds to an external magnetic field in 
physics and to a curvature in context of the discrete geometry of graphs. 

Let us return to the Bose-Einstein condensation. Let I be a positive integer 
and we consider T; of l|1.21fl and the finite graph obtained by the periodic 
boundary condition. Consider the adjacency matrices A^ and A r p. Let v be 

the periodic potential of Theorem 11.81 v is a periodic potential for the finite 

(v) 

periodic system on T^' as well. Set 

hi = - A rW + v, h = E-A r +v, (3.9) 

where E^ (resp. E) is the supremum of the spectrum of A r ( P ) — v (resp. Ar — v). 

1 

Now we consider the particle density. Set 



Pi( z ) = 



1 / ze- f3h 



r (p) 
1 / 



tr 



l2 ( T l P ') \ l-ze-P h 
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Due to translational invariance, 

1 f 1 / ze-Wto \ 

p(z) = lim pi(z) = — — 3 / — - tr Z 2 (Vo) dp. (3.10) 

i^oo (27r) d J[_ 7riOT] <i |r | \l - ze-l 311 ^ J 

As there exist constants Cg and Cl^such that 

c (i) 1 < < c (2)i 

^ x + 1 — z ~ 1 — ze - ' 311 — ^ x 

for positive x we have the following lemma. 

Lemma 3.2 Let -E(p) be the largest eigenvalue of the matrix Ar(p) — v(p). 
(i) Suppose the following integral is finite. 

' dp <oo. (3.11) 



[ — 7T)7r] u E-E(p) 

Then p(l) is finite and for any z n and l n satisfying lim n z n = 1 and lim n i n = oo . 

Urn Pl Jz n ) = p(l) (3.12) 

n — >oo 

(z?) Suppose that the following limit is infinite. 

lim / — ; dp^oo. (3.13) 

z/ij^^ E-E(p) + l-z y V ; 

Then, for any large positive p, there exists z ( < z < 1) such that p = p([3). 

Note that for any periodic potential v 

v{p) = v\ Vo 

Thus by abuse of notation we identify v{p) and v. 

v(p) — v 

H3.13J) of the above Lemma implies the absence of Bose-Einstein condensation 
while the case (i) of the above Lemma suggests Bose-Einstein condensation for 
the mean particle density p greater than p(l). We show that Ij3.13|l holds if the 
dimension v of our periodic lattice is one or two and that (|3.11|) is hnite if if the 
dimension v is greater than or equal to three. 

Lemma 3.3 

E = E(0)>E(p). (3.14) 
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The proof is same as periodic Schrodinger operators on R^. For the detail of 
proof, see Chapter XIII, 16 of |T^] . 

Lemma 3.4 In a neighborhood of p = 0, 

\E-E(p)\ <C\\pf = cj2pl (3-15) 

fe=i 



Proof. We set 

Re(i r (p)) = ~ (l r (p) + (A r (-p) 

When the quasi- momentum p is sufficiently small, Ke(Ar{p)) and hence, Ke(Ar(p))- 
v are Perron-Frobenius positive matrices and we have a positive vector as a 
unique eigenvector for the largest eigenvalue. Then, 

E(p) > sup (/, (Re(2 r (p)) - v)f) P{Vo) = sup (|/| , {Re(A r (p)) - v} |/|) P(Vo) 
ll/ll=i ll/ll=i 

(3.16) 

where |/| is a vector 1 2 (Tq) with the component \fj\ (j E Tq). 

The matrix element of Re(Ar{p)) is the real part of (|3.4() and (|3.5(l . 

[Re(Mp))h = cos^-(p)[Re(2r(0))]«. 

When the quasi-momentum p is sufficiently small, we have a small positive 
constant C such that 

cos^-Cp) > (1-CbH 2 ). 

Thus 

sup (|/|,{Re(2r(p))-4|/|)p (Vo ) 

ll/ll=i 

> (1 - C \\pf) sup (|/| , {(A r (0)) - v} \f\) P(Vo) 
ll/ll=i 

= (l-C\\p\\ 2 )E(0). (3.17) 
This inequality implies Ij3.15|l . End of Proof. 

Lemma 13.41 shows the divergence of the integral (|3.13|l if v is one or two. 
Next we consider the case v > 3. To show Proposition 13.61 below, we use a 
graph analogue of the diamagnetic inequality. 
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Lemma 3.5 For any function f onVo, 



r f ^-p iE -A(p)+v)A <{\f\,e-^ E - A ^+^\f\) . (3.18) 

V / y V / Vo 



Proof. Recall that the abosulute value of matrix elements of A(p) is same as 
that of A(Q) By the Trotter-Kato product formula, we have 



-l3(E-A(p)+v) 



' V N^o 



e N 



N 



N 



Vo 



< Jim |/| , 



e « 



V AT y 



TV 



l/l 



Vo 



(|/|, e -^-*°>+«> |/| 



Vo 



(3.19) 



End of Proof. 



Proposition 3.6 If E — E(po) for somepo ^ 0, there exists a diagonal unitary 
W on 1 2 {Vq) such that 



WA(p)W* = A(p - p Q ) 



(3.20) 



in a neighborhood of po 



Proof. Suppose that E = E(0) = E(po) and / is the unit eigenvector for the 
largest eigenvalue of A(p) — v. 

(A(p)-v)f = Ef, ||/|| = 1. 

By the diamagnetic inequality, 

-P(E-A( )+v) |jA < x (3 21) 

/ Vo v ' V 

Thus by differentiating (|3.21|) . 

(1(0)-«)|/|=^|/|. 

It turns out that |/| is the Perron Frobenius eigenvector of A(0) — v and all the 
components of |/| are positive. Now we define S a ( a € To ) and a diagonal 
unitary V via the following equations: 



l=(/,e-^ s -^)+")/) Vo <(|/|,, 



f a = e a '\f a \, W = diag(e- <a -). 
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By definition,/ = W* |/| , WvW* = v and 

(WA(p )W* -v)\f\=E\f\. (3.22) 

We claim that 

WA(p)W* = A(p - p Q ). (3.23) 

First, due to l{3~8l) . 

[WA(p)W*]tf = e^^* 4 "^ [!(())];,• . (3.24) 
Taking the real part of Ij3.22|) . we have 

(B(p )-v)\f\=E\f\, (3.25) 

where 

B{p) = Rc{WA{ Pa )W*). 
The matrix elements of B(p) is given by 

[fl(p)]« = cos(%(p) + * - ^)[2((%. 
Consider the following inner product: 

(|/| , (B(p ) - v) |/|) ;2(Vq) = E (|/| , |/|) ja(Vo) = E (3.26) 

As |/| is the Perron Frobenius eigenvector of A(po) — v, we have two identities: 

^ (cos(%(p ) + Si - 6^1(0)}^ ~ %) \f\ t \f\ 3 = E, (3.27) 

£ ([M0)]a - l/l, \f\j = E (3.28) 

y 

As these two identities must be valid simultaneously and all the components of 
|/| are non- vanishing, 

cos(%(p ) + Si - 5j) = 1. 
This implies sm(8ij(po) + Si — 5j) = and 

e V=I(««(p)+*i-«j) = e V=T(0y(p-Po)+(e«Cpo)+«<-«;O) = e V=T(9«(p-po)) ( 3 _29) 

The equations (eqn:z761) and (|3.29() suggest the claim of Propositior ]3.6l 
End of Proof. 

The following is a corollary of Propositior QSO 
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Corollary 3.7 (i) The number of po satisfying E — E(po) is finite. 

(ii) The integral ]S.11\) is finite if and only if the following integral is finite in 

a neighborhood of p = 0: 



Remark 3.8 (i) For Schrodinger operators with periodic potential on Euclidean 
spaces, E — E(p) holds only at the origin p = See J^j. 

As A(p) — v can be interpreted as a Schrodinger operators with a magnetic 
field on the finite graph Tq, degenracy E — E(p) can happen even if p ^ 0. Such 
examples are presented explicitly by Yusuke Higuchi and Tomoyuki Shirai in J^j. 
(ii) The above Provosition V3.b\ and Corollary \3. 7| are valid for any dimension v. 



Next we compare E—E(p) when v is the degree of the periodic lattice (i.e. when 
h = d— A = — Ar) and other v. Now E d , E d (p) stand for E and E(p) when v is 
the degree d and E v , E v (p) for those of other periodic potential v. By definition, 
E d = 0. We also denote h(v) = E - A r + v and hi(v) = E^ - A (p) + v . 

Lemma 3.9 Let Q = fl(i)(i S V) and QS l > — be the positive periodic 

ground states for h(v) andhi{v) satisfying the following normalization condition: 

{n,n) lHTa) = m)\ 2 = |v | , (rt l \rt% HVo) = |Vo] 

Then, 

(fn,h(v)fn) P{v) = \f(i)-fU)\ 2 nmU) (3.3i) 

for any f in l 2 (V) and 

(/n w ,fti(«)/n w )p( Vl )= E \f(i)-fU)\ 2 ^ l) m {l) (.i) (3.32) 

for any f in 1 2 (Vq). 

If f satisfies the twisted boundary condition with quasi- momentum p, we have 

(/n, (e - A T {p) + «)/n) P(Vb) = J2 |/( J ) - eiM?) /(i)f mm (3-33) 

where h^(v), andS^ are obtained by Aj? — v andSo with the periodic boundary 
condition. 9ij(p) is defined in \3. b)) when i and j are connected by an edge 
not belonging to £q and when i and j are connected by an edge in £q, we set 

OiM = o. 
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Proof. To derive 1)3. 33[) , divide (|3.32(1 by the volume of T/ and take the limit 
of I to the infinity. Thus we have only to show the identity (|3.31|1 which can be 
obtained by direct calculation as follows, (c.f. jH].) 

First we consider the case of d = v. In this case ft(i) = 1 and by definition, 

(fn,(d-A r )fQ) P{v) = e (\m\ 2 +\fU)\ 2 -7U)m-f(i)m) 

(i,j)ef 

E i/(0-/(i)i 2 - ( 3 - 34 ) 

Next we set w = v - d. Then ,-(d - A r )fi = (E + w)Q. Using $T5fy . 
(jn,(E + d-A r + w)fCl) P(y) 

= (fn, (d - A r )fn) lHV) - (/n, f(d - A r )n) nv) 

E \f(iW)-mm\ 2 

- E {i/wi 2 (i°(oi 2 - "(one?)) + \f(f)\ 2 (mtf - mm))} 

(i,i)e£ 

E (l/WI 2 + l/(i)l 2 -(7«/(i) + 7(i)/«)}W(i) 

(».j')e£ 

E \f(i)-m\ 2 mm- (3.35) 

End of Proof. 

Proof of Theoren il.81 By Lemma f3. 91 . we have 

E v (0)-E v (p)= inf [ ' J > " , , 

Eievo 1/(01 Mi)) 2 

E (l7lF£ <p) 1/(0 - /(j)l 2 

E d (0) - E d (p) = -E d (p) = mf -J^Mo __ . (3. 36 ) 



Set 



M 



^ Eiev i/wr 

supi n(i) 



infj ri(i) 
Due to (j3~3fi|) , we obtain 

M- 1 (E d {Q) - E d (p)) < E v {0) - E v (p) < M{E d {0) - E d (p)). (3.37) 

This inequality shows that finiteness of the integral i|3.3U[l in a neighborhood of 
p = for the discrete Laplacian is equivalent to that for the Schrodinger operator 
h with a periodic potential v. This completes our proof of Theorerr ll.81 End 
of Proof. 
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